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ABSTRACT 


This  papor  considers  the  equations  relating  errors  in  inertial 
guidance  system  measurement  instruments  to  errors  in  missile  burnout 
position  and  velocity  and  consequently  to  target  misses.  The  icIaHon- 
ship  is  expressed  in  terms  of  three  linear-  differential  equations 
with  trajectory-dependent  coefficients  and  forcing  functions  which 
depend  on  the  measurement  instrument  errors.  The  forcing  functions 
are  presented  in  forms  sufficiently  general  to  accomodate  errors 
that  might  arise  in  widely  different  types  of  instruments,  both 
accelerometers  and  gyroscopes.  Furthermore,  the  forcing  functions 
can  be  evaluated  for  arbitrary  instrument  orientations.  Programmed 
on  a  digital  computer,  these  equations  allow  swift  computation  of 
target  misses  for  given  component  tolerances  and  component  orienta¬ 
tion  and  consequently  facilitate  the  optimization  of  instrument 
orientation.  Such  a  program  has  been  used  extensively  by  the  authors 
for  gyro  orientation  schemes  and  general  error  studies  for  a  number 
of  ballistic  missile  inertial  guidance  systems. 
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A  DIGITAL  COMPUTER  PROGRAM  FOR  A  GENERALIZED 
INERTIAL  GUIDANCE  SYSTEM  ERROR  ANALYSIS 


1.  INTRODUCTION 

The  purpose  of  this  report  is  to  describe  and  explain  the  use  of  a  digital 
computer  program  which  is,  at  present,  on  the  1103A.  This  program  relates  inertial 
guidance  system  component  errors  to  errors  in  position  and  velocity  at  burnout  and 
consequently  to  target  misses.  The  program  has  two  prime  uses,  one  being  to  perform 
complete  guidance  system  error  analyses  and  the  '~>ther  being  to  study  the  optimiza¬ 
tion  of  instrument  orientations.  In  the  past,  error  analyses  and  design  studies 
for  inertial  guidsuice  systems  have  been  performed  either  by  hand  or  on  an  analogue 
computer  and  for  specific  instrument  (accelerometer  and  gyros)  configurations. 

Aside  from  the  obvious  need  for  speed  and  permanency,  the  need  was  felt  for  a  pro¬ 
gram  in  which  the  forms  of  the  error  functions  arising  from  the  instrument  errors 
were  sufficiently  general  to  accomodate  many  different  instrument  configurations 
(particularly  gyros)  where  the  instruments  are  arbitrarily  oriented  with  respect 
to  a  stable  platform.  It  is  felt  that  this  program  is  sufficiently  general  to 
satisfy  this  need.  The  authors  wish  to  extend  their  thanks  to  Messrs.  Frank  Meek 
and  Vilas  Henderson  for  their  very  capable  efforts  in  setting  the  problem  up  on  the 
computer. 

The  inputs  to  the  program  are  the  missile  trajectory  and  the  corresponding 
miss  coefficients;  and  the  inertial  component  errors  such  as  accelerometer  errors, 
gyro  drifts,  platform  misalignments  as  well  as  the  appropriate  angles  describing 
the  orientations  of  the  instruments.  The  outputs  are  the  position  and  velocity 
errors  at  burnout  and  the  corresponding  target  misses.  Section  5  concerned  with 
a  description  of  the  use  of  the  program,  particularly  the  input  of  instrument  errors. 
Sections  3  and  4  give  the  derivations  of  the  equations  relating  burnout  position 
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and  velocity  errors  to  component  errors  and  the  derivation  of  the  forms  of  the 
forcing  functions,  representing  .he  instrument  errors.  Before  deriving  the  error 
equations,  the  coordinate  systems  used  will  be  discussed  in  Section  ?. 


Z.  COOK PI  NATE  SYSTEMS 

There  are  three  basic  coordinate  systems  used  in  this  problem.  Two  of  them 
are  represented  in  the  diagram  below:  the  XYZ  system  which  is  an  earth-centered, 
inertial  coordinate  system  (ECl)  and  the  xyz  system  which  is  also  inertial,  but 
centered  at  the  missile  launch  point  at  the  time  of  launch.  The  XYZ  system  is 
suen  that  the  Z  axis  is  the  earth's  spin  axis  and  the  X  axis  is  such  that  the  X-Z 
p.^ane  contains  the  launch  meridian  at  the  time  of  launch.  The  Y  axis  completes 
the  orthogonal  set.  In  the  xyz  system,  the  z  axis  is  the  local  vertical  at  the 
launch  point  at  the  time  of  launch,  the  x  axis  is  3own  range  in  the  direction  of 
the  target,  and  the  y  axis  is  the  lateral  axis  nor.ual  to  the  x-z  plane  and  completes 
the  right-handed  system.  The  angle  A  is  the  firing  azimuth  measured  southward  from 
the  launch  meridian  to  the  x  axis.  ^  and  ^  are,  respectively,  the  geodetic  and 
geocentric  latitudes  of  the  launch  point.  The  quantities  a  and  b  are,  respectively, 
the  semi -major  and  semi -minor  axes  of  tne  earth,  e  the  eccentricity,  and  H  the 
distance  from  the  center  of  the  earth  to  the  launch  point. 
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The  coordinate  transformation  relating  the  XYZ  and  the  xyz  systems  is 


X 

X 

y 

»  M 

Y 

z 

Z 

+  C 


(2.1) 


where  M  is  the  rotation  matrix. 


-cos  A  sin 

sin  A 

COB  A  COS 

M  « 

-sin  A  sin 

-cos  A 

sin  A  cos  ^  * 

(2.2) 

COS  fSj* 

0 

sin 

The  translation  vector  C  is 


C  =  R 


6in  6  cos  A 
sin  6  sin  A 
-cos  6 


( :;*3) 


The  third  coordinate  system  is  the  accelerometer  coordinate  system 
The  axes  of  this  system  are  the  axes  along  which  the  range,  lateral,  and  pitch 
accelerometers'  sensitive  axes  are  pointed.  The  accelerometer  coordinate  axes 
are  obtained  by  rotating  first  about  the  z  axis  through  an  angle  which  defines 
an  intermediate  system  =  z-  The  axes  are  then  rotated  about  y^  through 

an  angle  which  defines  an  intermediate  system  x^y^  =  ^ ±>z ?’  The  fl-nal  rota¬ 
tion  is  through  the  angle  £  about  the  x0  =  x.  axis. 
z--^z. 


At Jerome ter  Coordinate  System  Axes 
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The  coordinate  transformation  relating  the  accelerometer  coordinate  system  XA’  yA'  ZA 
to  the  xyz  system  is 


A l^roui'* tf'r  Or! rntatl  on 


koo647 
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A’  ,yA'  A 


An  optional  racanB  of 
the  rotations  are  defined 


rotating  the  accelerometers  has  been  provided  in  which 
as  in  Figure  2.04. 


(2.4; 


s  may 
rms 


First  rotation:  ©2  about  y  axis. 


FIGURE  2.04 


Optional  Accelerometer  Orientation 


Note  that  the  first  rotation  has  been  defined  0O.  The  matrix  of  rotations  is 


cos 

G1 

COS  ©2 

-sin  9^ 

-cos 

G1 

sin 

e2 

z  = 

sin 

G1 

cos  02 

cos  0^ 

-sin 

G1 

sin 

*2 

sin 

92 

0 

COD 

*2 

The  first  rotation  has  been  defined  to  be  ©2  because  by  setting  ^  =  0^,  -  ©2, 

i.  -  0  in  Q  (equation  2-5),  and  taking  the  transpose  of  Q,  namely  Q',  one  finds 

3 

Z  =  Q' .  The  computer  has  been  programmed  to  take  advantage  of  this  fact. 
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The  orientation  of  the  gyro6copeB  with  respect  to  the  accelerometer  and  plat¬ 
form  coordimte  systems  will  be  d:  scuBsed  in  Section  4,  which  is  concerned  with 
the  forcing  functions  of  the  error  equations. 


3 .  ERROR  EQUATIONS 

In  the  XIZ  system  let  represent  the  vector  of  missile  Inertial  accelerations 


r  1 


and  let  A^a  represent  the  vector  of  sensed  accelerat:  ons 


Futher,  let  denote  the  missile  position  vector 

(3-3) 

and  let  G  denote  the  gravity  vector 

X 

Y  (3-*0 
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where  gQ  is  gravitational  acceleration  at  the  launch  site,  Rq  is  the  distance  from 
the  center  of  the  earth  to  the  launch  site,  and  where  R  is  the  distance  from  the 
center  of  the  earth  to  the  missile.  Thus, 

R2  =  X2  +  Y2  +  z2  =  Pj*  ?1  (3.5) 

where  Pj'  denotes  the  transpose  of  the  vector  P^. 


Z 


FIGURE  3.01  M  ssile  Trajectory 


Then,  from  the  relationship 


AIa  +  0 


(3.6) 


one  obtains 
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Mj  •  MIa  ,  tc 


where  the  errors  in  inertial  accelerations  are 


*1  = 


is2 


and  where  the  errors  in  sensed  accelerations  sure 


AA 


la 


^a_ 


The  quantity  the  error  in  the  gravitational  vector,  is 


R  R 

*  ■  (3«o  -r^>  Fi  -  (s0  ±  ^ 


(3- 


(3-C 


(3-S 


(3.1 


wh.£rs 


Air 

AZ 


(3-1 


Equation  (3*10)  may  be  written  as 


AG  ■=  (36.  )  (RAR)  Pj  -  (g0  )  APj 


'°  R5 


(3.1 
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the  expressions  in  brackets  all  being  scalars.  The  scalar  may  be  written  as 

3-Y) 


R/*=XdX+Y£Y+Z£Z=  P 
Substituting  (3*13)  into(3»12)  yields 

,-.■8)  R  2  R  2 

*  ■  (3«o  4-Pi'^i)  Px  -  U0 

JK  n 

Substituting  ( 3 . 14 )  into  (3*7)  yields  the  error  equations  in  vector  fom 

R  2  R  2 

4  (go  3“  '  '  ^3so  3"  PI  ”  ^la 

3-9)  B  R 

This  equation  written  in  component  form  is 


(3.13) 


(3*lU) 


(3.15) 


3-10) 


3.11) 


^  +  e. 


3S, 


R  “ 

X(X£X  +  Y£T  +  ZdZ) 

V? 


*®Xa 


+  k 


Y(XAX  +  Y£Y 
R^ 


^  +  6, 


3s, 


Z(XdX  +  YdY  +  ZAZ)  =  AaZa 


(3-lo) 


or  three  simultaneous  linear  Becond-order  equations  which  may  be  solved  to  yield 
the  relationship  between  velocity  and  position  errors  (dX,£Y,&&)  and  (£X,AY,AZ) 
3.I;4)  and  the  errors  in  sensed  accelerations.  The  coefficients  in  the  differential 
equation  are  trajectory-dependent,  i.e.,  R,X,Y  and  Z  depend  on  the'  trajectory. 
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As  will  be  seen  in  the  next  section  on  forcing  functions,  the  errors  Aa^a,  Aa^a 
Aa^  are  also  trajectory* dependent  for  the  most  part. 

<uQ 

In  the  actual  computer  program  the  equations  are  solved  in  the  earth-centered 
inertial  coordinate  system  XYZ.  The  position  and  velocity  errors  resulting  from 
the  solution  are  printed  out  in  the  XYZ  system,  the  xyz  or  platform  system,  and 
the  yA, zA  or  accelerometer  system.  The  target  misses  are  also  printed  out,  but 
of  course  these  are  independent  of  the  coordinate  system.  In  actuality,  the  misses 
are  computed  in  the  XYZ  system,  the  miss  coefficients  having  been  appropriately 
transformed  from  the  xyz  system  to  the  XYZ  system. 


Since  the  sensed  accelerations  are  measured  in  the  x^, y^, system,  it  is  of 
interest  to  point  out  the  relationship  between  the  sensed  acceleration  errors 
Mja  in  the  XYZ  system  and  the  sensed  acceleration  error  AA^  in  the  accelerometer 
system.  It  is,  by  virtue  of  relations  (2.1)  and (2.4), 


(3-3 


from  which  it  follows  that 


Mja  -  (3-3 

since  the  inverse  of  QM,viz.  (QM)"''~  =  (QM)'  =  M'Q'  where  the  prime  denotes  transpose. 
Thus,  equation  (3*15)  becomes 


R  2  R  2 

+  (go  -f-  )  ^  -  <3g0  Pj'tfj)  Pj  =■ 

K  rt 


(3-: 


Equation  (3*19)  then  expresses  the  relationship  between  the  position  and  velocity 
errors  in  the  XYZ  system  and  errors  in  sensed  acceleration  in  the  accelerometer 
coordinate  system.  The  next  section  is  concerned  with  the  form  of  AA^  and  its 
relation  to  the  trajectory  being  used. 
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4.  FORCING  FUNCTIONS 

The  term  forcing  function  will  be  applied  to  the  vector  and  denoted 


by  the  vector  F.  Thus, 


F  =  M'Q'M 


(4.1) 


The  errors  In  sensed  acceleration  in  the  platform  system,  may  be  expressed 


MaA  -  h  -  *AaA 


(4.2) 


3-17) 


where  E^  is  the  vector  of  errors  arising  due  to  accelerometer  inaccuracies  such  as 
bias  and  scale  factor.  The  term  $A^  represents  acceleration  errors  due  to  angular 
errors  in  the  platform.  §  represents  the  platform  angle  matrix,  and  expresses 

the  spurious  accelerations  sensed  by  the  accelerometers  due  to  the  platform  mis- 
orientation  resulting  from  either  initial  platform  misalignment  or  from  gyro  drift. 
The  $  matrix  is 


3-10) 


^y  ~K 


C'.3) 


3.19) 


where  the  angles  , />  and  represent  rotations  about  the  x.,y*  and  zA  axes, 
x  y  z  n 

respectively. 


From  (2.l)  and  (2.4)  it  is  seen  that  A^  may  be  expressed  as 


A  .  =  QMAt 
aA  la 


(4.4) 


but  from  (3«6)  +  G  so  that 


\a  ■  W  (Aj  -  0) 


(4-5) 
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Thus,  from  (4.1)  and  (4.2)  above,  the  forcing  function  F  may  be  expressed  as 

F  =  M'Q’Ea  +  M’Q'$QM  (A  -  G)  (4.6) 

Section  5  is  concerned  with  the  general  forms  of  E^  and  $  due  to  accelerometer 
errors,  platform  angles,  and  gyro  drifts.  Before  proceeding  to  Section  5,  it  will 
be  helpful  to  describe  the  gyro  and  accelerometer  axes. 

Two  two-degree -of -freedom  gyros  have  been  assumed  and  their  orientation  restric¬ 
ted  to  cause  the  three  input  axes  to  form  on  orthogonal  set.  The  axes  h^, f^,  f,,  8X6 
the  axes  of  the  number  1  gyre  with  h^  the  input  axis.  The  axes  h2,h^, f^  refer  to 
the  number  2  gyro  with  h^  and  h^  the  input  axes.  With  respect  to  the  x^, y^, 
system  and  consequently  with  respect  to  the  xyz  system,  the  gyros  may  be  arbitrarily 
oriented,  subject  only  to  the  constraint  as  mentioned  above,  so  that  the  three  input 
axes  form  an  orthogonal  set.  Figures  4.01  and  4.02  illustrate  these  orientations. 


GYRO  NUMBER  1 
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GYRO  NUMBER  2 


*A 


ingle 

about 

when 

♦l 

z ^  axis 

*2  *  P2  *  0 

♦2 

h^  axis 

p2  =  0 

^2 

f ^  axi 3 

Gyro  Orientation 


Thus,  the  orthogonal  set  h^, is  related  to  the  ZA  system  by  the  rotation 

matiix. 


cos  *2  cos 

cos  ^2  sin 

sin  t2 

! 

( -cos  p  sin 

1  ^ 

(cos  P2  CCS  ♦1 

sin  P2  cos  ♦ 

R  = 

-sin  P2  sin  cos  ♦1) 

-sin  P2  sin  *2  sin  +  ) 

j 

( sin  P2  sin  ^ 

-cos  P2  sin  t2  cos  ^1) 

(-sin  P2  cos 

-cos  P2  sin  y2  sin  ♦1) 

cos  P2  cos  *2 
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FIGURE  4.03 


Gyro  Orientation 


Thus,  rotations  about  the  A  axes,  are  related  to  rotations  due 

to  gyro  drifts  about  the  h^h^h  axes,  say  ty  the  relationship 


|X1 

_ / 

i 

y 

=  R' 

*z 

L  zj 

*3_ 

where  R'  is  the  transpose  of  (4.7)*  It  can  be  shown  that  as  a  result  of  (4.8) 
the  $  matrix,  (4-3),  may  be  written  as 
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Thus,  the  forcing  function  (4.6)  may  be  written  as 
F  -  M’Q'Ea  +  M,Q'R'$hRQM(AI  -  0) 


]  where  j>h  is  the  matrix  .in  (4.9),  viz. 


(4.10) 


(4.11) 


It  should  be  pointed  out  that  the  gyro  arrangement  need  not  necessarily  be  two 
two - degree- of ^ -freedom  gyros.  Any  arrangement  wherein  there  are  drifts  about  three 
orthogonal  gyro  axes  may  be  used.  As  a  matter  of  fact,  many  cases  of  non-orthogonal 
axes  may  be  treated  by  using  one  set  of  values  of  the  angles  in  the  matrix  R  for 
the  number  1  gyro  and  in  an  independent  computer  run  use  another  set  of  angles  in 
R  for  the  number  2  gyro. 


Certain  types  of  accelerometers  and  velocity  me-. era  are  such  that  the  plane 
in  which  the  pendulous  element  moves  and  which  contains  the  sensitive  axis  may  be 
at  any  arbitrary  angle  with  respect  to  the  x^,  y A,  system.  Denote  the  axes  of 

the  meter  itself  as  1,2,3,  where  1  is  me  sensitive  axis  and  axes  1  and  2  define 


the  plane  of  motion  of  the  pendulous  element. 


Axis  3  completes  the  orthogonal  set 


and  is  the  axis  about  which  the  pendulous  element  rotates 


2 
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Figure  4.05  illustrates  two  possible  orientations  of  the  2  and  3  axes  for 
the  pitch  accelerometer. 


FIGURE  4.05  Accelerometer  Alignment 


Due  to  the  physical  nature  of  the  meters  and  the  unequal  accelerations  along 
the  and  axes,  there  may  be  some  preferred  orientation  which  minimizes 

the  cross-axis  acceleration  errors.  For  this  reason,  provisions  in  the  program 
have  been  made  for  arbitrary  orientations  of  the  meters  about  the  sensitive  axes. 
The  cross-axis  accelerations  (accelerations  along  axes  2  and  3)  are 


- 

a„ 

0  cos  0  sin  0 

a  . 

2x 

X  X 

xA 

a_ 

_ 

-sin  9  0  -cos  8 

a  . 
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<< 

V! 

yA 

a. 

-sin  Q  cos  0  0 

a  A 

Hz. 

z  z 

—  — 

zA 
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where  a2x  is  the  acceleration  along  axis  2  of  the  x  or  range  accelerometer,  a^„ 
the  acceleration  along  axis  2  of  the  y  or  lateral  accelerometer,  and  where  ©^  1b 
a  rotation  about  the  sensitive  axis  of  the  range  accelerometer,  ©  is  a  rotation 

y 

about  the  sensitive  axis  of  the  lateral  accelerometer  and  ©  is  a  rotation  about 

z 

the  sensitive  axis  of  the  pitch  accelerometer.  The  angles  are  described  in 
Figure  4.06.  Although  the  accelerometers  used  here  for  illustration  have  a  pendu¬ 
lous  element  moving  in  the  1-2  plane,  the  relationships  given  in  (4.12)  and  (4.13) 
may  be  used  for  other  types  of  instruments  for  the  study  of  cross  accelerations. 


/ 

i 


FIGURE  4.06 


Accelerometer  Alignment 
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The  next  section  describes  the  forms  of  the  elements  of  E^  and  $  and  4^  in 
terms  of  missile  accelerations  and  actual  component  tolerances  and  the  manner 
in  which  they  are  identified  and  entered  into  the  program. 


5.  SIMULATION  OF  ERROR  SOURCES 

The  purpost  of  an  inertial  measurement  unit  is  to  provide,  as  its  output,  a 
true  indication  of  any  acceleration  (excepting  gravity)  the  unit  is  experiencing. 

To  accomplish  this,  the  function  relating  accelerometer  output  to  input  should 
be  well  known;  and  the  orientation  of  each  accelerometer,  with  respect  to  some 
coordinate  system,  should  be  well  defined  both  in  space  and  time.  When  these 
conditions  are  not  satisfied  an  acceleration  error  appears,  and  it  is  the  purpose 
of  the  veetor 

F  =  M*Q'Ea  +  M'Q'R'^RQMUj  -  G)  ( 

to  evaluate  this  error. 

Broadly  speaking,  acceleration  errors  can  be  catalogued  into  two  classes.  The 
errors  arise  because: 

CLASS  1:  The  accelerometer  output-input  relation  is  not  well 

known  (e.g.,  scale  factor  uncertainty)  or  is  non¬ 
linear,  or  the  accelerometer  is  misaligned  with  respect 
to  the  platform, 

CLASS  2:  The  stabilized  platform  on  which  the  accelerometers 

rest  is  initially  misaligned  or  has  rotated  from  its 
correct  orientation. 

The  two  classes  can  be  associated  with  the  two  parts  of  the  vector  F  as  was 

discussed  in  Section  4.  Thus,  to  simulate  ClasB  1  errors  must  be  evaluated; 

while  for  Class  2  errors  must  be  formed.  The  rest  of  this  section  will  be  con- 

h 

cerned  with  the  simulation  of  each  class  with  attention  being  focused  first  on 
Class  1  errors. 
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ClaaB  1  ErrorB 

Error  sources  which  fall  in  this  category  are: 


(a)  Accelerometer  misalignments  with  respect  to  the  platform. 

(b)  Uncertainty  in  zero  setting  (null  shift  or  bias). 

(c)  Linear  scale  factor  uncertainty. 

( d)  Non-linear  terms  in  the  input-output  relation  -  proportional  to  square, 
cube,  cross-products  and  absolute  value  of  acceleration. 


To  simulate  Class  1  errors  the  following  general  vector  has  been  programmed 
for  evaluating  E^. 


(5-2) 


where  the  are  3  x  3  matrices, 
functions  of  acceleration. 


ard  the 


are  vectors  who Be  elements  are 
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(5.3) 
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rv 
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T? 
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a2x  * 

a  . 
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layA  1  • 

a  k 

yA 
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Equation  (5-3)  continued 


EQ0  E8l  E82 

had'  °xa 

E90  E9i  E92 

hx1  •  %A 

E83  E84  % 

la2y |  •  syA 

+ 

E93  E94  E95 

la3y|  '  ay* 

_E86  E87  Eai 

l&2z 1 '  azA 

i 

E96  E97  E98_ 

ia3z|‘V 

r_ 

n 

r 

r_ 

_  ~1 

r  n 

E100  E101  E102 

E110  Elll  E112 

+ 

e103  eio4  e105 

^10 

+ 

E113  E114  E115 

^11 

_E106  E107  E108 

_^116  E117  Ell8_ 

axA’ ayA* azA  exe  range,  lateral,  and  pitch  accelerometer  sensed  accelerations;  and 

a2x’ a2y’ a2z’ a3x' a3y' a3z  ^  defined  in  equations  (4.12)  and  (4.13)  with  a^, ayA, 
a^  again  defined  to  be  sensed  accelerations.  Note  that  provision  has  been  made 

undefined. 

Thus,  regardless  of  the  error  source,  if  it  can  be  expressed  as  a  function  of 
acceleration  it  can  be  simulated  in  this  program. 

By  use  of  this  general,  vector  the  error  sources  listed  above  can  be  simulated 
in  the  following  ways. 

a.  Accelerometer  Misalignments  with  Respect  bo  the  Platform 

Figure  5*01  shows  the  rai so  dentations  of  the  accelerometers  that  are  to 

be  simulated.  3  is  the  angular  rotation,  in  radians,  of  the  x  accelerometer 
xy 

toward  the  y.  axis,  0  is  the  rotation  of  the  z  accelerometer  toward  the  x. 

A  2X  A 

axis,  and  similarly  for  the  others. 


for  simulating  any  function  of  acceleration  by  leaving  Ag^Ag^, 
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FIGURE  5.01 


Accelerometer  Misalignments 


The  acceleration  errors  resulting  from  such  misorientations  are  simulated  by  the 
1  allowing . 


0  6 

Kxy 

axA 

EA  “  ^E2^E2  = 

P  0 

yx 

<V* 

V 
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b.  Uncertainty  In  Sero  Setting  (Null  Shift  or  Blag) 

If  a  ,a .a„  represent  the  shift  or  bins  in  the  zero  settings  of  the 

X  J  *4 

three  accelerometers,  then  the  error  in  the  indicated  output  can  be  represented 
by 


ea  -  *1  - 


c.  Linear  Scale  Factor  Uncertainty 


(5-5) 


Uncertainties  in  scale  factor  setting  give  rise  to  errors  which  can  be 
simulated  by  setting  the  elements  of  to 


€x  0  0 
°  €y  0 

0  0c 

z 

then 

*A  *  ^E2^E2 


(5.6) 


where  «  ,c  ,c  are  the  uncertainties  in  the  range,  lateral,  and  pitch  accel 
x  y  z 

erometers,  respectively. 


d.  Non-linear  Terms  in  the  Input-Output  Relation 

To  illustrate  the  simulation  of  non-linearities,  the  hypothetical  relation 


aout 


a3 

‘la  A 


(5-7) 


will  be  assumed  for  each  aceelerometer.  Then 


TR-59"0000-0u6)i7 
Page  23 


represents  the  error  in  the  accelerometer  output.  Such  errors  can  be  simula- 
ed  by  setting 


„  .th 

wnere  ^jLi’^Si  8X6  i:ne  ior  t-n,;  1  accererometer .  moie  exampj.es 

of  such  simulations  will  be  presented  in  Section  6. 

In  all  of  these  illustrations,  in  fact  whenever  the  vector  is  being  used, 
it  is  necessary  to  specify  the  angles  5^,  ly  or  ^ePenciin6  on  the  option 

being  used  to  specify  accelerometer  orientation  (  see  Section  2 ) . 


Class  2  Errors 

Error  sources  considered  to  fall  in  this  category  are: 


(a)  Initial  platform  misalignments. 

(b)  Fiatform  rotations  due  to  platform  elasticity. 

(c)  Platform  rotations  due  to  constant  or  acceleration-dependent  gyro 
drift  rates. 


T?-y>.oooo-ou<34Y 

Page 


When  pln'fonu  rotations  are  considered  to  be  small,  the  acceleration  error 
con  be  found  from 

M  =  -  ?  x  A  (5.1C 

where  ^  denotes  the  angular  rotation  vector  of  the  stabilized  platform  with 
respeat  to  its  proper  orientation.  Expanding  the  vector  cross-product  in  a 
coordinate  system  defined  by  the  axes  h1 , h0, h,  (see  Figure  4.03),  and  putting 
the  result  into  matrix  form,  yields 

*S,  *  Vh  -  *hMA  -  -  °>  <5-U 

where is  defined  by  equation  (4.1l)  and  R,  defined  by  equation  (4.7),  rotates 
accelerations  along  accelerometer  coordinates  into  the  h^, hg,h  axes.  Multiplying 
by  R'  rotates  the  acceleration  errors  into  accelerometer  coordinates,  and  multi¬ 
plying  by  M'Q'  completes  the  rotation  into  ECI  coordinates  in  which  the  differential 
equations  are  solved.  To  simulate  Class  2  errors  the  rotation  vector,  must  be 
evaluated.  To  accomplish  this,  a  general  vector  G  has  been  programmed  from  which 
p  is  evaluated  by  either  of  two  options. 
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P-10) 


}  ^  MGiAGi 
i  =  2 


(5-iM 


in  which  the  are  3x3  matrices  and  the  A 
are  functions  of  acceleration. 


are  vectors  whose  elements 
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As  in  the  vector  E^,  provision  has  been  made  for  simulating  any  function  of 
acceleration  desired  by  leaving  AG6  and  Aq-,  undefined.  The  acceleration  vectors 
appearing  in  G  are  defined  by 

M3) 


mn  mr AAifl.  t-f 
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The  matrix  S  is  identical  to  R  except  that  P2  must  "be  replaced  by  ^  (see  equation 
4.7).  Reference  to  Figures  4.01  and  4 02  showB  that  the  effect  of  the  matrices  S 
and  R  is  to  rotate  accelerometer  accelerations  into  the  axes  defining  gyros  1  and 
2,  respectively.  Note  that  the  matrix  R  is  also  used  in  the  vector  F  to  obtain 
accelerations  along  the  three  gyro  input  axes,  h^,  h  ,  h^;  or,  from  another  viewpoint 
it  is  used,  as  shown  in  equation  (4.9),  to  obtain 

$  =  R'$hP  (5.1-7) 


where  $  is  the  matrix  containing  rotations  about  platform  x,y,  z  axes.  It  is  impor¬ 


tant  to  keep  in  mind  the  two  roles  that  the  matrix  R  has  when  simulating  ui&ss  2 
errors.  This  will  be  demonstrated  in  the  following  examples. 


Initial  Platform  Mi  sails 


Normally,  initial  platform  misalignments  are  given  about  platform  x, y, z 

axes  as  <5  ,1  By  setting  t,  =  =  0  in  R,  in  which  case  R  become, 

xo  yo  zo  i.  Ct  Ct 

the  identity  matrix,  the  vector  becomes,  from  equation  (4.8), 


(5-13) 


and  =  $h  from  equation  (5-17)  ■  As  can  be  seen  from  figures  4.01  and  4.02, 
such  a  choice  of  angles  does  indeed  cause  the  axes  to  lie  along  the  x, y,z 
axes.  The  simulation  is  completed  by  choosing  Option  1,  i.e.,  /  =  G,  and 


S' 


[ITU 
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0  ' 
v-> 

0 

_ 1 

—  — " 

^XO 

°u 

E 

6 

ryo 

C12 

rzo 

1 

J 

—  — 

in  equation  (5.I5). 


I 

(5-19) 


If,  however,  the  initial  platform  mis&ligniaants  were  given  about  axes 
other  than  x,  y,  z,  say  accelerometer  axes,  which  were  rotated  through  angleB 
$1»  Jg,  L,  then  the  simulation  would  be  effected  by  setting 


/*•  — 

r  “1 

«i 

Xo 

^10 

*2 

- 

*2 

,  *  G,  and 

°n 

m 

^20 

_*3_ 

.V 

._^30 

(5-20) 


In  effect,  the  h^  axes  have  been  placed  parallel  to  the  accelerometer  axes, 
acceleration  errors  computed,  and  then  rotated  back  to  platform  x,  y,  z  axes 
and  finally,  as  is  done  in  all  cases,  into  15CI  coordinates. 


b .  Platform  Rotations  Due  to  Platform  Elasticity 

If  the  platform  rotates  when  accelerated,  the  indicated  acceleration 
errors  can  be  found  by  choosing  the  angles  and  such  that  the  h^ 

axes  lie  along  the  axes  about  which  the  platform  rotates;  and  using  the 
vector  G  to  simulate  the  expressions  which  relate  rotations  about  the  axes 
to  accelerations  along  these  axes.  For  example,  if 


^1  “  Klahl  +  K2Qh3  +  K3ah2ah3 


h 


=  0 


(5-21) 
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the  simulation  vould  be  effected  by  choosing  Option  1,  i.e.,  ^  =  G,  where 


“K1  °  K2_ 

af3 

"0  k3  o' 

af  3  ’  ^2 

°  =  MG2AG2  +  *04*04  “ 

0  0  0 

V 

0  0  0 

®h2  *  ®h3 

0  ^  K5 

L.  J 

*h3 

0  0  0 

®h3  ’  af3 

NOTE:  a^1  -  (see  Figures  4.01  and  4.02) 


In  the  two  examples  given  thus  far,  the  matrix  R  haB  served  as  a  means  of 
obtaining  accelerations  along  the  h^  axes  with  no  mention  of  gyro  axeB.  The  h, 
axes  have  been  used  as  an  intermediate  set  which  in  some  cases  are  more  convenient 
to  use  than  platform  x,y, z  or  accelerometer  axes.  In  the  case  of  gyro  drift  rates, 
the  h^  axes  will  be  used  to  represent  the  three  axes  stabilized  by  the  gyros.  The 
matrix  R,  though  still  used  in  the  vector  F  to  obtain  accelerations  along  the  hj 
axes,  is  looked  upon  more  as  a  means  of  obtaining  accelerations  along  the  axes  of 
gyro  2.  The  matrix  S  is,  of  course,  used  to  obtain  accelerations  along  the  axes 
of  gyro  1.  Because  the  simulation  of  gyro  drift  rates  is  dependent  on  specific 
gyro  orientations  and  drift  rate  equations,  examples  will  not  be  given  until 
Section  6  in  which  a  specific  gyro  system  will  be  simulated. 

Before  proceeding  to  Section  6,  it  should  be  pointed  out  that  the  vectors 
and  G  can  be  used  simultaneously,  but  that  only  one  optic l  may  be  used  for  evalua¬ 
ting  p  for  any  given  simulation. 

6.  DETAILS  OF  COMPUTER  PROGRAM 

To  use  the  program  it  is  necessary  to  provide  a  tape  on  which  are  recorded  the 
sensed  accelerations,  position,  and  range  from  the  center  of  the  earth  of  the  missili 
as  a  function  of  time,  all  in  ECI  coordinates.  It  is  advantageous,  but  not  neces¬ 
sary,  to  have  recorded,  at  burnout,  the  missile  velocity  components  in  ECI  coordin¬ 
ates.  Once  such  a  tape  is  made  available,  the  input  required  for  evaluation  of 
error  sources  is: 

( a)  General  Input 

Azimuth  angle  A 


1. 
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(«0 

General  Input  (continued) 

i«. 

Gravity  components  in  ECI  coordinates  if  (E^q,  1 

°r 

are  to  be  computed  (see  examples  for  use  of  this  feature) 

(  b) 

For 

Vector  E 

A 

1. 

Angles  ^  or  ei.©2 

2. 

Wea  “  (VVa)  or  (a3x’a3y’ a32'  a” 

being  used 

3- 

Elements  of  the  matrices 

(c) 

For 

Vector  G 

1. 

Angles  \|r?,  P1,P2 

2. 

Whether  Option  1  or  2  is  being  used 

3- 

Elements  of  the  1^,  matrices 

10*  11'  12  ‘ 


Elements  and  angles  not  specified  are  assumed  to  be  zero.  All  angles  should  he 
given  as  positive  numbers  in  degrees, 

Figure  o  *  01  she  vs  iiow  ou't.pu't  ini' onus,  on  is  pi*£  ^  Blocks  1  s.nd.  2  reproduce 

whatever  _nput  has  been  provided  to  the  program,  while  Blocks  3  and  ^  present  the 
errors  computed  as  a  result  of  this  input.  Position  and  velocity  errors  are  printed 
in  accelerometer,  platform  x,  y,  z,  and  ECI  coordinates;  and  if  burnout  velocity  has 
been  recorded  on  tape,  perturbed  position  and  velocity  vectors  are  printed  in 
ECI  coordinates  which  can  be  used  in  free-f light  programs  to  determine  the  target 
miss.  The  target  miss,  in  terms  of  down-range  and  cross-range  miss,  is  computed 
in  Block  h  if  mise  coefficients  are  provided. 

The  term  "coupled  or  uncoupled"  is  used  to  denote  the  form  of  che  differential 

equations  being  solved.  That  is,  one  has  the  option  of  solving  equation  (3>l6)  with 

g  1  0  or  g  =  0,  in  which  cases  the  equations  are  said  to  be  coupler  or  uncoupled. 

0  0 

Uncoupling  the  equations  eliminates  whatever  effect  £G  (equation  3*12,  gravity 
computation  errors)  may  have  on  the  final  position  and  velocity  errors. 


ECI 

Mi sg  Coefficients 

«D 

M„ 


Trajectory 

Identification 

Number 


Errors  in 

Accelerometer 

Coordinates 


Errors  in  Launch 
Point  Coordinates 
(x,y,z) 

Errors  in  ECI 
Coordinates 


Elements  of  p 
Elements  of  '/ 

ECI  Position  and 
Acceleration 


ECI  Perturbed 
Position  i  'Aslocity 
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gQ  (ft/ sec2) 

R-,  u  (ft) 
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Though  not  indicated  in  the  output  format,  provision  has  been  made  in  the 
program  whereby  initial  conditions  can  be  placed  on  &.,&f  and  £2. 

Figure  602  indicates  the  manner  in  which  the  computer  input  is  utilized  and 
how  the  program  operates.  A  test  for  non-zero  matrices  is  made  thereby  eliminating 
time  consuming  operations  with  matrices  that  are  identically  zero.  ^te  the  follow 
ing  about  the  marker  in  which  output  information  is  printed: 

(a)  Block  1  is  printed  only  once  for  any  given  trajectory. 

(b)  Block  2  is  printed  only  once  for  any  given  case. 

(c)  Block  3  only  is  printed  at  optional  printout  times. 

(d)  Block  4  is  printed  only  at  burnout. 


7.  EXAMPLES 

To  demonstrate  the  use  of  this  program,  several  examples  of  error  sources  that 
have  been  simulated  with  this  program  will  be  presented. 

Class  1  Errors 

Won-Linear  Terms  in  the  Input  -  Output  Relation 

The  accelerometers  used  in  one  guidance  system  are  of  the  pendulous  type 
shown  In  Figure  4.04  and  are  oriented  on  the  platform  as  shown  in  Figure  4.06, 

Three  axes,  1,2,3  are  used  to  describe  the  accelerometer,  with  axis  1  the  sensitive 
axis.  The  equation 

aout  -  Co  +  »1  +  Clal  4  C2  Kl  al  +  C3  N  al  +  C3a2al  +  cl,  N  al 

+  c5»3  ♦  -  (cA  ♦  CjlgJSi  +  c3  |e2)  Si  +  C's2gl  <  cu|e3]  S;L 

-•  *  C5«i2) 


(7.1) 
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hae  been  used  to  relate  output  to  Input  acceleration  for  any  given  accelerometer 
where  a^  ■  aA,  the  acceleration  along  the  sensitive  axis,  and  g^,  g^,  are  gravity 


components  appearing  along  the  three  axes  at  launch, 
accelerometer  output  is  then 


The  error  in  the  indicated 


**  *  aout  '  al  (T,2) 

For  demonstration  purposes,  and  because  the  rest  of  the  terms  are  simulated  in  a 
similar  manner,  only  terms  involving  C^,C^  will  be  simulated.  For  these 

terms 


EA  ”  ^1  +  +  ^7^7  + 
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En 
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0 

0 

C3*_ 
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azA_ 
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L. 

0  ck 

a2x  axA 
a2y  ayA 
_a2z  azA„ 


ce  0  0 

5x 
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p 

0  C_  0 

5y 

•3. 

yA 

0  0 

a3. 

L  5zJ 

L  zA_ 

where 


"E10 

C3x  lg2x)  ®xA  +  C’3xg2xgxA  +  C5x®xA 

E11 

3  - 

C3y  +  C3yg2ysyA  T  C5ygyA 

E12 

Cjz  ^g2zl  gzA  +  ^3zg2zgzA  +  ^5igzA 

(7.3) 


(7.U) 


in  addition,  the  angles  or  roust  be  specified  to  define  the  accel- 

'•r .z.-t'-r  f.enultlve  axes  with  respect  to  the  platform  (x,y,z)  axes; 

it.  :  t.}.»  angler.  9  .ft  ,9  must  be  given  to  define  axis  2  of  each  accelerometer  with 

X  /  Z  0 
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respect  to  the  (x^, y^, zA)  system.  If,  for  example,  axis  2  were  along  the  zA  axis 
for  the  range  accelerometer  and  along  the  -yA  axis  for  the  pitch  accelerometer 
(see  Figure  1».06),  one  would  set  6 


90°,  6, 


180 


Except  for  special  orientations,  the  gravity  components  g^,  gg,  g^  are  complica¬ 
ted  functions  of  the  angles  specified  above,  thus  necessitating  a  lengthy  computa¬ 
tion  to  evaluate  equation  ( 7-M  •  To  eliminate  this  computation  (especially 
desirable  for  optimization  studies  where  the  angles  are  being  varied)  the  computer 
has  been  programmed  to  perform  this  operation  at  t  =  0.  The  computation  is  made 
as  follows: 


r 

/ 

1 

r  2  n 

r  31 

~\ 

E10 

1 

gxA 

®xA 

«xA 

Eu 

-  -  < 

^2 

V 

+  ME3 

2 

®yA 

+  mei* 

g3 

gyA 

+  •  •  *  > 

_®zAJ 

2 

l®^AJ 

j 

In  other  words,  the  acceleration  vectors  appearing  in  EA  are  computed  at  t  =  0 
in  terms  of  gravity,  matrix  multiplications  and  summations  are  carried  out,  and 
the  negative  of  this  sum  put  into  the  cells  used  for  the  elements  ^11' ^12' 

To  use  this  feature,  gravity  components  existing,  at  launch,  and  expressed  in  ECI 
coordinates,  must  be  provided. 


Figures  7>01  and  7-02  are  examples  of  the  output  obtained  when  C  ,  C'  and 

are  simulated  individually.  Thi  elements  E,^  and  E.,„  have  been  machine 
yz  JO  12 

computed  by  use  of  the  feature  described  above.  The  following  input  was  required. 


Figure  7*01  - 
Figure  7.02  - 
Figure  7.02  - 


3  x 
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x  10'6 

9x  ■  20° 

=  315° 

£2  = 

3y 
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x  10'6 

ey  "  10° 

ix  =  315° 

t2  =  ^5° 

5z 

Eu8  -  0.3 

x  10'6 

=  3150 

*2  -  1*5° 

■  n,  gravity  conjionents  and  miss  coefficients  were  provided. 
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In  another  guidance  system  the  accelerometers  are  oriented  such  that  their 
sensitive  axes  lie  along  the  platform  (x,y,  z)  axes,  i.e.,  *  tg  «  *  0*  The 
equations  used  to  describe  the  error  in  their  outputs  are 


An  .  "  K  +  K.  a  +  K0  a  2  4  K„  a  3  +  K.  la  a  I 
xA  ox  lx  x  2x  x  3x  x  4x  !  x  xl 

*  V  +  Vy  +  V/  +  Kiyay3 

'  Koz  *  KiA  *  *  k3z^3  +  KUz  IV.I 


(7.6) 


Since  a  a  =  a  a  for  this  system,  the  equations  can  be  simulated  as  follows: 
X  Z  X  z 


EA  =  “ei  +  ^2*32  +  ME3AE3  +  ME4AE4  +  ”eA7 
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I  (7.7) 

'where  ©  =  90°,  ©  =  270°. 

x  z 
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aa 
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Class  2  Errors 

Platform  Rotations  Due  to  Platform  Elasticity 

The  platform  in  one  guidance  system  is  such  that  when  subjected  to  accelerations 
it  tends  to  rotate,  with  the  rotation  being  directly  proportional  to  acceleration. 

The  rotations  are  about  platform  (x,y, i)  axes  where 
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From  equation  (4.0)^  when  ik  *» 

*  (+ 


Po 
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6 

(7.9) 


With  such  a  choice  of  angles,  the  error  source  la  simulated  by  choosing  Option  1, 
i  .e.,  j)  =  G;  and 

ro 

1 

G  =  ”02*02  =  '  ° 


The  only  input  required  to  simulate  this  error  source  is 


9  -ax 

0 

0 


0 

0 


a  1 
1'3 

*U2 

>3J 


(a 


*3 


ahl  =  ax> 


(7.10) 


f>  =  U, 


\1  r f 

22  x ' 


‘26 


Figure  7-03  is  an  example  of  ouch  a  simulation  where  a  =  a  =  0.2  x  10  .  Note 

that  since  all  angles  are  zero,  none  have  been  printed. 


uy  ru  ul  i.  I  i/  net  uc  & 

A  typical  gyro  can  be  defined  by  three  orthogonal  axes,  S,  IW,  9W,  as  shown 
in  Figure  7.04. 


S 
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The  acceleration-dependent,  drift  ratte  about  bxcb  IW,  OW  can  be  represented  for  a 
certain  system  as 


^IW  =  Bsaiw  +  °saiw  "  "  KIWaSaIW 


*^OW  =  BSaOW  *  ■  KOWasaW 


(7.11 


B. ,0. ,BG.  are  ball,  global,  and  ball-giobal  unbalances  along  the  i  axis,  and 
Ki  =  (Kg  "  •*,s  ^e  difference  in  compliance  between  the  spin  and  i  axis. 

A  platform  can  be  stabilized  by  two  Buch  gyros  oriented  with  respect  to  the 
platform  as  shown  in  Figure  7.05. 


GYBO  1 


7R  -  so  -  (X>CO  - 006-  7 
Fage  -1 


cKsncnly  used  to  define  the  (-S)  axis  of  ©  to  1  (pitch  gyro), 
lie  sieved;  thus,  three  drift  rate  equations  suffice  to  describe 


'^:wi  =  3si®:*i  "  ssiaiwi  '  3omas:  "  Kiwiasiaiwi 

^IW2  c  4  GS2£IW2  "  EGIW2aS2  *  KIW2aS2aTW2 

^0W2  =  BS2a0W2  +  BOW2aS2  *  KOW2aS2&OW2 


This  orientation  can  be  simulated  in  the  error  analysis  prograa  by  choosing 
Vp  and  defining  the  gyro  axes  as  shewn  in  Figure  ~t . 06  (refer  to  Figures 

h .01  end  h.02). 


GYRO  1 


GYRO  2 
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It  is  necessary  to  add  the  2?0°  offset  in  defining  p^  because  in  the  simulation  ^ 
is  used  to  define  axis  f^  with  respect  to  the  (-x)  axis,  where  p^  is  measured 
counterclockwise  about  h^.  When  defined  in  this  manner,  the  angle  p,  defined  in 
Figure  7*05*  is  equal  to  7.  Equations  (7. 12)now  become 

^1  “  ^lVL  +  °81ahl  "  ^IWlaf 2  '  KIWlaf2Bhl 

*^2  “  BS2ah2  +  GS2fth2  "  ^IWS^S  "  (’ 

*^3  "  BS2°h3  +  BOW2af3  ’  KOW2af3ah3 


Under  the  assumption  that  a  2  o,  then  a. 

■Jf 


f3  "  Vl 


0,  and  equations  (7.13)  reduce  to 


^1  “  BGIWlaf2  (dropping  the  minus  sign) 

^2  “  BS2ah2  +  °S2ah2 
’U  —  11  - 

^3  "C2“h3 

t 

Equations  (7*13)  can  be  sit alated  by  choosing  Option  2,  l.e.,  ^  =  J'  Gdt  and 

o 

G  *  *^2*02  +  >\l3AG3  +  *1c4AG4  +  *^05 


0 

0 

— \ 

0  i 

(BS1  +  °81^ 

0 

_BGIW1 

-B0IW2 

^GS2  +  °82^  ° 

A  + 

*G2 

0 

0 

0 

_B0W2 

0 

BS2 

0 

0 

0 

0 

1 

0  0 

■ 

0 

0 

< 

Su 

0  0 

AG4  + 

0  0 

0 

*05 

1 

0 

0  ~KCM2 

0  0 

0 

(1 


(7 
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I  Figure  7.07  chows  the  errors  that  ar^se  due  to  where  0 ^  -  0.3014  x  10"^, 

I  0^  »  322°  corresponding  to  p  »  52°. 


|  Constant  Gyro  Drift  Rates 

1 

1  Constant  gyro  drifts  are  simulated  by  choosing  the  same  angles  and  option  and 

i 

|  setting 


(7.13) 


G  “  " 


r .  n 

^10 

^20 

jSo 


(7-14) 


Though  drifts  due  to  gravity  are  not  subtracted  out  at  launch  in  this  Gystem, 
provision  har  been  made  to  evaluate  the  elements  Gio,Gu'Gi2  ln  terros  oi  gravity 
components  Just  as  the  elements  E^,E^,E^  can  be  evaluated  in  the  EA  vector. 


(7-15) 


FIGURE  7-07  Computer  Output  for  BGIW1 
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